Mon. Not. R. Astron. Soc. 000, [TJ{8] (0000) Printed 22 November 2012 (MN WF&t style file v2.2) 



The influence of star formation history on the gravitational 
wave signal from close double degenerates in the thin disc 



Accepted . Received 



<n: 

O- Shenghua Yu and C. Simon Jeffery 

^SJ 1 Armagh Observatory, College Hill, Armagh BT61 9DG, N. Ireland 

o 

(N 
< 

o 

43 
6 



ABSTRACT 

The expected gravitational wave (GW) signal due to double degenerates (DDs) 
in the thin Galactic disc is calculated using a Monte Carlo simulation. The number 
of young close DDs that will contribute observable discrete signals in the frequency 
range 1.58 — 15.8 mHz is estimated by comparison with the sensitivity of proposed GW 
observatories. The present-day DD population is examined as a function of Galactic 
star-formation history alone. It is shown that the frequency distribution, in particular, 
is a sensitive function of the Galactic star formation history and could be used to 
measure the time since the last major star-formation epoch. 

Key words: stars: white dwarfs — stars: evolution — stars: binaries: close — stars: 
formation — Galaxy: evolution — Galaxy: structure 
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1 INTRODUCTION 

Theoretical studies have indicated that double degener- 
ate (DD) stars with very close orbital separation should 
be the most abundant and prominent sources of gravita- 
tional w ave (GW) radiation at frequencies between 0.1 and 
30 mHz JPeters fc Mathews! 1 19651 : lLandau fc Lifshit zl 1 19751 ; 
lEvans et all 1 19871 ). In addition to other binaries, they will 
form a GW confusion background due to the large num- 
ber of sources per fr equency resolution bin l|Hils et al.lll990l : 
iHils fc BendeiTl2000h . 

In general, the number of DDs per equal-interval fre- 
quency bin at high frequency in any Galactic popula- 
tion drops as a function of frequency. For a frequency 
resolution better than ~ 3 x 10 -8 Hz and frequencies 
>1 mHz (i.e. orbital periods ^ 33min), a GW experi- 
ment may detect only one or no individual sources per 
resolution bin, ma k ing the source poten t ially resolvabljj 
jEvans et all Il987); IHils fc Bender] bood: INelemans et af 
l200ll : iRuiter et al.H2009l : lYu fc Jeffery||2010l ). Resolved DDs 
identified from GW detections would allow us to study the 
relation between t heir properties and t he dynamic nature of 
globular clusters (|Willems et al.l [20071 ). the nature and po - 
tential outcome of the DD merging process l|Webbinklll984 ), 



and their association with the formation and evolution of a 
galaxy. 

Since DDs in a galaxy are the outcome of stellar for- 
mation and evolution, a key question is how the star forma- 
tion (SF) history, the distribution of initial orbital param- 
eters and the evolution of their main-sequence progenitors 
influence the DD population in terms of final distribution 
with orbital period (GW frequency), separation and chirp 
masfl This question can be addressed using the method 
of binary-star population synthesis (BSPS) to investigate, 
for example, the SF histo ry alone. In our previous papers 
|Yu fc Jeffervl |2010| . I20TH ), the influence of star formation 
and evolution on the DDs in a galaxy was characterised by 
the use of a parameteric SF rate and synthetic SF contribu- 
tion functions. 

The SF history in a galaxy is likely to be more com- 
plicated than approximations used in previous studies, for 
example, a si ngle star-burst, c onsta nt SF, or (quasi-) ex - 
ponential SF (INelemans et all l|200ll ): IRuiter et all (|2009h ; 
lYu fc Jeffervl (|2010l )). Observations agree that the SF rate 
in the Galactic disc is likely to b e declining but is episod - 
ically enhanced (|Maiewskilll993l ; iRocha-Pinto et ai]|2000h . 
The aim of this paper is to show how the GW signal from 



1 If there is only one DD in a resolution bin and the GW signal 
from this DD is greater than a given detection threshold (e.g. 
S/N=3 in this paper), we define the DD to be a resolved source. 
If the GW signal from this DD is less than the same detection 
threshold (e.g. S/N=3 in this paper), we define the DD to be a 
potentially resolved source. 



2 Chirp mass (M) of a binary is defined as Ai = 

( mi+m 2 ) ^ m2 ) 2//5 wnere m i an d rn-2 is the mass of each 

component of the binary respectively. A general relation between 
the GW frequency / and the orbital period P or b in the nth har- 
monic is / = p" . In particular, if a binary is in a circular orbit 

(eccentricity e=0), we have / = p 2 . 
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the DD population in the thin disc in the Galaxy depends 
on the SF history whilst keeping other factors in the model 
constant. 

We compute the GW signal from a simulated DD pop- 
ulation in the thin disc of the Galaxy in §[2] and compare 
the results for three different SF models in §[3] We compare 
our results with previous studies in §|4] and draw significant 
conclusions in §[5] 



2 THE GW SIGNAL FROM RESOLVED 
DOUBLE DEGENERATES IN THE THIN 
DISC 

We simulate the DD population in the thin disc of the 
Galaxy by adopting a simplified binary-star evolution and 
population-synthesis model (BSPS). Inputs include a star 
formation rate (SFR), an initial mass function (IMF), a mass 
ratio distribution P q , an orbital separation distribution P a , 
an eccentricity distribution P e , a me tallicity Z and a rep re- 
sentation of the thin disc structure l)Yu fc Jeffervl l20ich . If 
the orbital period of a binary is sufficiently long (i.e. >100 
yrs), the binary is assumed to behave as two single stars 
and will not form a DD. Inputs also include binary-star evo- 
lution physics, including mass-loss r ates, common - envelo pe 
ejection efficiency, etc., described bv lYu fc Jeffervl l)201(f ). 

To provide a benchmark calculation, the SFR is approx- 
imated by a quasi-exponential function 



S(t) = 7.92e" t/T + 0.09t M Q yr" 



(1) 



where t is the age of the disc and r = 9 Gyr, yield- 
an SFR w 3.5Mnyr _1 and a total disc mass (in 



mg 



stars) M d ~ 5.2 x 10 M© at the current epoch (t = 10 
Gyr). The value fo r the current SFR is consistent with re- 
ce nt observation s (Il5"iehl et al.l I2006T I. We adopt the IMF 
of iKroupa et all (1 19931), a constant P q ilMazeh et al.lll992l ; 
Goldberg fc Mazehlll994 ). and P a following iHanl lll998h ~and 



Griffin! l| 19851 ). P„ = 2e follows iNelemans et al.l (|200ll ). We 



adopt Z = 0.02 (solar) for all thin disc binaries. 

The simulation is started with a sample of 10 7 primor- 
dial MS+MS binaries, with total mass » 1.05 x 10 7 M Q . Ini- 
tial parameters for each binary are allocated on the above 
distributions using a Monte Carlo procedure. This proce- 
dure yields a total of ~ 4.93 x 10 4 DDs over a period of 
15 Gyr. We compute the birth rate, merger rate and total 
number, and record their individual properties (masses and 
orbital parameters) as a function of the age of the thin disc 
l|Yu fc Jeffervll201lh . 

Subsequently, a second Monte Carlo procedure is used 
to generate the orbital properties for a complete sample of 
Udd thin-disc DDs by interpolation on the initial DD sam- 
ple described above. The initial DD sample obtained from 
computations of stell ar evolution and th e complete sam- 
ple are ident i cal w ith lYu fc Jeffervl ()201ll ) , and similar to 
lYu fc Jeffervl l|2010h (the difference to the total number DDs 
in the thin disc is ~ 8%) . We identify four types of DD pair 
by the chemical composition of their degenerate cores, i.e. 
He+He, CO+He, CO+CO, and ONeMg+X where X = He, 
CO, or ONeMg. Since we use a simplified model, the crite- 
ria for identifying the WD c ore compositions are given by 
the initial-final mass relation (Han 19981: iHurlev et all 2000. 
12002ft . 



lYu fc Jeffervl (|201 ll ) give n d d for the four different types 
of DD in this simulation. These ndd DDs are distributed in 
a thin disc model where the stellar mass density is described 



Pd(R,z) 



Md 



4nh R h z 



- R/hR sech 2 (-z/h z ) M pc" 



(2) 



and where R and z are the natural cylindrical coordinates of 
the axisymmetric disc, and Kr = 2.5 kpc and h x = 0.352 kpc 
are the radial scale length and sca le height of the thin disc 
i|Sackettlll997l ; iPhleps et afll200d ). This structure and the 
thin disc mass (in st ars) are consis t ent w ith t he dynamical 
param eters given by iKlypin et ail (|2002h and iRobin et al.l 
i|2003h . 

The DDs are sorted to obtain the number distribution 
by orbital frequency, to calculate the total strain amplitude 
^2 h 2 from all DDs in each frequency bin, and hence to ob- 
tain the GW strain amplitude-frequency relation due to the 
thin disc DDs. We adopt a frequency resolution element 
Af = 3.17 x 10~ 8 Hz, corresponding to one year of observa- 
tion with a GW detector. 

Figure [T] illustrates the influence from different SF 
epochs on the predicted GW signal. The expected sensitiv- 
ities for S/N=3 for the proposed GW detectors eLISA and 
LAGRANGE ( or eLISA-type mission) in a one year m ission 
are also shown ((Larson et al.l[200ol : IConklin et al.ll201lh . 

From this figure, we identify a turning point, or critical 
frequency / c , around 0.001 Hz. For frequencies / < / c , the 
amplitude of the GW varies approximately inversely with 
frequency, forming a 'continuum'. For frequencies above the 
critical frequency f c , the amplitude becomes a discrete sig- 
nal. 

The existence of the critical frequency is strongly asso- 
ciated with the product of the number distribution of DDs 
with respect to the frequency and the size of the frequency 
element, i.e. ^fi-Af. When / sC /«., i.e. ^f-Af > 1 , 
there must be one or more DDs per resolution element Af, 
implying that n(f) can be represented by a continuous func- 



tion. When / > f c , i.e. 



dn(f) 
df 



Af < 1 , some frequency bins 



will not contain any DD. Rather, there will be one DD in 
i = ( dn ^P A/) -1 frequency bins, i.e. in one frequency bin 
between / and / + iAf. 

These results are c onsistent with an analytic approach 
fe.g. lEvans et all < |l987h ). whereby 



fc = 2.977 



M 

Mr,, 



-5/11 



yi- 



3/11 



Af 

10~ 8 



3/11 



mHz, 



(3) 

in which the chirp mass M of a binary is assumed to be 
an independent variable and v = J f(f)df is the total birth 
rate of DDs. 

Figure[T] shows that the continuum of the thin disc GW 
signal consists mainly of signal from early SF (i.e. 0.0 — 6.0 
Gyr). Late SF (i.e. 8.0—9.95 Gyr) contributes most of the 
discrete signal in the high frequency band and only a small 
fraction of the continuum. Sources with log / > —2.3 princi- 
pally comprise CO+He and He+He DDs formed during the 
last epoch of SF. In the range —3 < log / < —2.3, different 
types of DD will be mixed, with CO+CO and ONeMg+X 
DDs dominating the high amplitude signal and CO+He and 
He+He DDs contributing to the low amplitude. We are pri- 
marily interested in these resolved sources since their strain 
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Figure 1. The contribution of the current DD populations from different star-formation epochs to the gravitational- wave amplitude 
(log h) in the quasi-exponential star-formation model. Solid (cyan) and broken (magenta) lines represent the eLISA and LAGRANGE 
sensitivity with S/N=3 for a one year mission, respectively. 



amplitudes exceed the sensitivities of the proposed GW de- 
tectors. 

In addition, distance strongly affects the GW signal, 
both continuum and discrete. The structure of the thin disc 
implies that the distance to the majority of DDs will be 
large, but there will be a few nearby (local) DDs. Thus, in 
our simulations, there are few isolated strong signals due to 
DDs at d < lkpc (e.g. Fig.QJb): log / = -2.88 and log / = 
-2.98). 



3 COMPARISON OF DIFFERENT STAR 
FORMATION MODELS 

We have discussed the GW signal with a 'Quasi- Exponential 
SF' rate. In order to see the influence on the GW signal 
from different SF models, we have calculated the GW sig- 
nal assuming two further models for the SF rate. These are 
'Instantaneous SF', in which a single star burst takes place 
only at the formation of the thin disc with a constant SF 
rate of 132. 9M yr" 1 from to 391 Myr, and 'Constant SF', 
in which SF occurs at a constant rate of 5.2Mq yr" 1 from 
to 10 Gyr. All three models produce a thin-disc star mass 
of ~ 5.2 x 10 10 M Q at the thin-disc age t = 10 Gyr. 

The differences between the signals in these models arise 
from DDs with different progenitor properties. According to 
binary-star evolution theory, main sequence stars with mass 



M ^ 2M© can leave low-mass white dwarfs (WDs) with 
core mass <0.5 M © if th ere is no significant mass accretion 
(e.g. iHurlev et ail (|2002l )). Stars with M « 2 - 10M Q can 
develop WDs with massive degenerate cores (>0.5 M©) and 
heavier elements (i.e. CO or ONeMg). Since the evolution 
time-scale of a star is a strong function of its initial mass, 
the birth rates of different WD types and the GW signal are 
sensitive to different SF epochs. For example, the signals at 
high frequency from He+He and CO+He DDs are sensitive 
to SF between 6 — 9 and 8 — 9.4 Gyr respectively (Fig.[T}. Our 
results show that a constant SF model produces a distinctly 
higher strain amplitude from He+He DDs (—22.2 < log/i < 
-23.0, -2.8 < log / < -2.0) than both other models. The 
quasi-exponential model produces less high-frequency GW 
radiation than the constant SF model. 

However, a very different discrete GW signal at / > 
0.001 Hz is produced by the instantaneous SF model. The 
bottom panel in Fig. [2] shows only a few GW signals from 
close DDs at / > 0.00158 Hz (log / > -2.8). Since the ma- 
jority of close D Ds merged in the pa st 10 Gyr (see birth and 
merger rates in lYu fc Jeffenl (|201lD ). only a few close DDs 
now remain. This implies that short-period DDs producing 
the discrete GW signal are sensitive to recent SF (see Fig.[T]). 

To represent possible future observations, the sensitiv- 
ities of the proposed eLISA and LAGRANGE missions are 
shown in each figure. In our present simulations, the number 
of (one-year) resolvable DDs are 11880, 22920 and 16240 for 
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Figure 2. Resolved and potentially resolved DDs in the different star-formation models. From top to bottom, the SF models are quasi- 
exponential, constant and instantaneous (§0. Solid (cyan) and broken (magenta) lines represent the cLISA and LAGRANGE sensitivity 
with S/N=3 for a one year mission, respectively. For clarity, the systems shown represent a small sample (one DD in twenty) of the 
complete thin-disc simulations. 



the quasi-exponential, constant, and instantaneous SF mod- 
els respectively, in which about 2690, 9020 and 80 could 
be resolved by LAGRANGE with S/N=3, and about 8010, 
19820, 3840 could be resolved by eLISA with S/N=3. These 
systems are shown in Fig. [21 The proportion of each type of 
DD is shown in Table[T] The proportions of resolvable DDs 
detectable by LAGRANGE are quite different from those 



given by the models because of the frequency response of 
the detector. LAGRANGE is more sensitive to high-/ DDs, 
of which the highest proportion are found in the constant 
SF model, which contains the highest fraction of DDs from 
recent star formation. 

Figure|3] shows the distribution of the GW strain am- 
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Table 1. Proportion of each type of resolved DD in each SF model. 





quasi-exp 


constant 


instantaneous 


types of DDs 








He+He 


71.0% 


78.9% 


74.6% 


CO+Hc 


22.6% 


16.2% 


17.1% 


CO+CO 


5.5% 


4.1% 


7.8% 


ONeMg+X 


0.9% 


0.8% 


0.4% 


total 


11880 


22920 


16240 


resolved by LAGRANGE with S/N=3 


2690 


9020 


80 


resolved by eLISA with S/N=3 


8010 


19820 


3840 



plitude and frequency from all resolvable DDs in each SF 
model. This figure can be understood analytically. 

The change of orbital separation a of a DD in a cir- 
cular orbit with respect to time due to GW radiation is 
i|Landau fc Lifshit z|[l975h 



64^3 _5 „ ,2 -3 

- — G c uM a 

5 



(4) 



where a is the orbital separation; G is the gravitational con- 
stant; c is the speed of light; fj, and M are the reduced mass 
and total mass of the DD, respectively. Combining Eq.0] 
with Kepler's third law P 2 rb /a 3 = 47r 2 /GM and the rela- 
tion between orbital period and GW frequency P or b =2//, 
we have the change of the GW frequency of a DD per unit 
time 



i 96 5/3 



c V /3 X 5/3 / 11/3 , 



(5) 



where M = ^ 3//5 M 2//5 is the chirp mass (see footnote 2). 
The inverse / _1 indicates the GW radiation timescale or 
lifetime. 

Assuming that the birth rate of the DDs at frequency 
/ is v(f), the number of DDs per unit frequency may be 
written as the product of birth rate and lifetime; using ap- 
propriate units 



Mf) 

d/ 



Kf) 



'-G 5 / 3 c~ 5 n s / 3 M 5 / 3 / n / 3 



^O.O546(A^/M )- 5/3 (/,(/)/yr- 1 ) (//Hz)" 11/3 Hz" 1 . 

(6) 

Hence the number distribution has contributions from the 
birth-rate v{f), which is important when significant num- 
bers of young DDs are being injected into the population at 
high frequencies, and from the orbital decay lifetime which 
gives the / _11//3 distribution and which is important when 
short-period DDs have all evolved from old long-period DDs. 

To illustrate, we assume that v(f) and M are inde- 
pendent of frequency /. If M = 0.5 Mq and v = 0.004 
yr _1 (average values estimated from simulations with quasi- 
exponential SF in the frequency range lO^-lO -1 ' 6 Hz), the 



number probability of DDs 



d«C0 

d/ 



0.0007/ 



-11/3 



(df = 1 



yr _1 ). This relation is plotted in Fig. [3] and has the same 
slope as the simulation with instantaneous SF (approxi- 
mately a (^-function) for log / > —2.8. 



The differences between the number distributions of re- 
solved DDs in the three SF models can be understood sim- 
ply in terms of the SF history. The DD birth rate u{f, t) can 
be expressed as v(f,t disc ) = / () tdisc C,(f,t diBC - t si )S{t B! )dt sf 
where tdisc is the age of the thin disc, t s t is the time of star 
formation, C, is the contribution function, and S is the star 
formation rate. C, provides information about the number 
and frequency distribution of DDs from a given epoch, and 
is independent of th e SF history (see §2.1 and Figs. 6-7 in 
lYu fc Jeffervl (|201ll ) for details). The DD birth rate is ob- 
tained by convolution with 5*. 

If S is a (5-function, the birth rate of DDs after 10 Gyr is 
a very weak function of /. Most present-day close DDs would 
be long-evolved DD systems, so their number distribution 
follows the power- law relation /~ 11/3 (Fig.|3]). However, if S 
is continuous over 10 Gyr, contributions at every SF epoch 
must be included and will increase the number of DDs at 
high frequency. 

Thus, assuming that (a) all other inputs to the BSPS 
model are valid (or at least that systematic errors can be 
eliminated by other means) and (b) there is no significant 
contamination by resolved sources from other locations, it 
is evident that GW observations of DDs could distinguish 
between different SF histories. In particular, the frequency 
above which the distribution develops a / _11 ' 3 gradient is 
indicative of the time since the last major SF epoch. 



4 DISCUSSION 

The impact of different SF histories provides one reason 
for the differences between calculations of DD number dis- 
tributions and GW strain am plitud e -frequency relation by 

lI |1987| ; INelemans et~ail l200ll ; 



different authors (Evans et al 



iRuiter et all [2009 



Yu fc Jeffervl |2010| ). especially where 



these authors adopt similar approximations for the physics 
of binary star evolution. The method of interpolation (see 
§[2} is important for predicting the absolute number of re- 
solved DDs but not crucial for the relative numbers in differ- 
ent models. In the present simulation, we only interpolate to 
produce a complete Galactic sample from the more limited 
BSPS sample. There is no extrapolation and the interpola- 
tion errors decrease the predicted GW strain and the number 
of r esolved DDs by no mor e than ~30%. 

INelemans et all l|2004h estimate that there may be w 
2.13 x 10 8 close detached DDs in the Galactic bulge and disc, 
with « 2.3 x 10 7 semi-detached DDs, assuming an age of 13.5 
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Figure 3. Histogram of the GW signal distribution from resolvable DDs, including resolved and potentially resolved DDs (top panels) and 
from resolved DDs by the two GW detectors (middle and bottom panels) with respect to strain amplitude (left ) and frequency (right ) in 
different SF models, with A log h = 0.1 and A log / = 0.05. Blue represents the instantaneous star formation (SF) model, red, constant SF, 
and green, quasi-exponential SF. The straight black line represents the number-frequency relation: dn ^j^ ~ 0.0007/ -11 / 3 , df = 1 yr _1 . 
Top panels: no correction for detector sensitivity, Middle: relative distributions assuming an eLISA-type detector with S/N=3, Bottom: 
relative distributions assuming a LAGRANGE-type detector with S/N=3. All number fractions are normalized to the total numbers of 
resolved DDs in each model (Table 1). 



Gyr and a SF history following iBoissier fc Prantzosl |l999). 
We find » 2.4 ± 0.5 x 10 8 DDs in the bulge and thin disc, 
with « 5 x 10 5 semi-detached DDs f w 0.3% of all DDs in the 
thin disc) llYu fc Jeffervll2010l.l201ll). Th e average birth rate 
of DDs given by iNelemans et alj l|2004l ) (1.58 x 10~ 2 yr -1 ) 



is slightly less than our result (2.4 x 10 yr ). However, 
we give a much lower average birth rate for semi-detached 
DDs « 5 x 10~ 5 yr"" 1 , compared with w 1.7 x 10 



-"yr" 1 



by INelemans et all (|2004l ). The local space density of semi- 
detached DDs in our quasi-exponential SF model of the thin 
disc is ~ 6.2 x 10 -7 pc -3 , slightly less than the constraint 
of 1 — 3 x 10~ 6 pc~ 3 obtai ned from the Sloan Digital Sky 
burvey l|Roelofs et al.ll2007T I on the assumption that semi- 
detached He+He and CO+He DDs look like AM Canum 
Venaticorum binaries. 
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iRuiter et all (|2010h estimate there to be » 2.48 x 10 7 
DDs with orbital periods in the range 200 s to 5.6 hr in the 
Galactic disc, including ~ 1.62 x 10 7 semi-detached DDs 
(~65% of all DDs). They assume a constant SF rate (4 
Mq yr~ ) in the disc for 10 Gyr. The average birth rate of 
DDs in this period range is ~ 2.48 x 10~ 3 yr~\ while the av- 
erage birth rate of semi-detached DDs is ~ 1.62 x 10 -3 yr _1 . 
The latter is significantly higher tha n our result. The lo - 
cal density of semi-detached DDs in IRuiter et all l|2010h . 
2.3 x 10 -5 pc -3 , is substantially higher than both our 
value and the observation. The formulae adopted to sim- 
ulate Galactic structure may contribute a significant frac- 
tion of the differences in estimates of the l o cal d ensity 
of DDs; iNelemans et all (|2004h . IRuiter et ail l|2010f ). and 
ourselves adopted density distributions p cc sech(z/h z ) 2 , 



p cc e 



(l-z|/M 



and p oc sech (— z/h z ) respectively, where 
p is mass density, z the height of the thin disc in the cylin- 
drical coordinates, and h z the scale height of the thin disc 

llYu fc Jefferj201Ch 

Nissanke et all <|2012T ) recently confirmed that thou- 



sands of detached DDs may be detected by space-borne GW 
detectors with 5 and 1 Mkm arms, but there may be only 
a few ten to a few hundred detections of interacting DD 
systems. This is attributed to an assessment of detection 
prospects, based on iterative identification and subtraction 
of bright sources with respect to both instrument and con- 
fusion noise. Our results also imply that detached DDs will 
dominate the GW signal from the total Galactic DD popu- 



lation in the GW frequency range of 10 



■10" 



' Hz. How- 



ever, they also argue that the SF history of the disc is an 
important factor and may affect the number ratio of the 
detached and semi-detached DDs detected by a GW obser- 
vatory. 

Since the GW background may be affected by sources 
including DDs beyond the Milky Way, an estimate of con- 
tamination from the Local Group would be of interest. 
Observations indicate that SF history in the Local Group 
is complicated. For instance, the Large Magellanic Cloud 
(LMC) most likely experienced a dramatic decrease of SF 
following an initial star burst some 12 Gyr ago. SF resumed 
some 5 Gyr ago has since proceeded with a low average rate 
of roughly 0.2 Af Q yr _1 , pea king roughly 2 Gyr, 0.5 Gyr, 
0.1 Gyr and 0.012 Gyr ago (|Harris fc Zaritskvl l2009l ) . The 
Small Magellanic Cloud (SMC) has a similar history to the 
LMC, but with a different look-back time for the initial star 
burst (~ 8.4 Gyr) an d with minor star bursts roughly 2.5, 
0.4, and 0.06 Gyr ago (|Harris fc Zaritsky||2004l ). Our results 
(e.g. Fig.QJ show that recent (2 Gyr) SF has a significant 
contribution to the GW signal from DDs in the Galactic thin 
disc. Assuming that DD contribution functions are the same 
for both the Galaxy and for the LMC and SMC, the GW 
signal from the thin-disc DDs may well be contaminated by 
DDs in other Local Group galaxies. This should be explored 
in the future. 



5 CONCLUSION 

Monte Carlo simulations have been carried out to synthesize 
the GW signal from DDs based on a binary-evolut ion model 
ijHanll 19981 : iHurlev et alJkool lYu fc Jeffervll2010h . different 
star-formation models and the present thin disc structure 



of the Galaxy. The GW signal consists of a continuum at 
low frequency and a discrete signal at high frequency. The 
continuum arises mainly from DDs with main-sequence pro- 
genitors formed in the early disc (i.e. between and 8 Gyr). 
Recent SF (i.e. between 8 and 10 Gyr) gives rise to a signifi- 
cant fraction of short-period DDs corresponding to the GW 
signal at high frequency, which is mostly resolved and which 
is thus sensitive to the recent star-formation history. At the 
highest frequencies, where the orbital-decay timescale gov- 
erns the evolution, we have shown that our Monte Carlo sim- 
ulations reproduce the distribution expected from analytic 
approximation. Both the intercept of the high-frequency dis- 
tribution of the resolved DD population and the frequency 
below which the distribution drops below that predicted by 
the orbital-decay timescale are sensitive to recent star for- 
mation. 

Since the high-frequency GW signal is dominated by 
DDs from recent star formation, future work should explore 
whether GW observations can be used to constrain the time 
since the most recent major burst of star formation in the 
Galaxy. 
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